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A  COMPUTER  ASSISTED  METHODOLOGY  FOR  ALIGNING 
THE  AN/ASN  43  FOR  DOPPLER  (AN/ASN  128) 
INSTALLATION  AND  CALIBRATION 


INTRODUCTION 

Uith  the  procurement  of  a  doppler  navigation  system  by  the  US  Army  for 
Installation  on  certain  Army  aircraft,  it  is  imperative  that  maximum 
fidelity  be  attained  from  the  aircraft  heading  reference  system.  This 
report  deals  only  with  the  alignment  of  an  AN/ASN  43;  however,  the 
methodology  described  is  appropriate  for  aligning  and  calibrating  other 
magnetic  reference  systems. 


ALIGNMENT 

A  detailed  method  for  aligning  the  magnetic  flux  system  is  contained 
in  the  AN/ASN-128  Doppler  Installation  and  Acceptance  Testing  Manual.1  This 
section  will  only  serve  to  supplement  and  enhance  the  base  document . 

In  initiating  an  alignment  procedure,  a  precise  flight  test  course 
should  be  established  for  use  throughout  the  entire  procedure.  Many  factors 
exist  which  must  be  considered  in  establishing  a  course.  Several  of  these 
factors  are  listed  below: 

1.  Select  check  points  (CPs)  that  can  be  easily  identified  at  a 
distance  of  several  kilometers  when  flying  the  test  course. 

2.  Select  a  CP  that  is  easily  identified  on  a  large  scale  map  to 
assure  the  greatest  accuracy  in  plotting  its  coordinates. 

3.  Design  the  test  course  so  that  the  CPs  selected  are  close  to 
magnetic  North-South  and  East-West  headings.  This  will  help  in  resolving 
errors  by  holding  the  number  of  possible  variables  to  a  minimum. 

4.  The  CPs  should  be  located  approximately  15  to  20  kilometers 
apart.  A  graphic  presentation  of  an  ideal  flight  test  course  is  depicted  in 
Figure  1. 


1  US  Army  Avionics  Research  &  Development  Activity.  Installation  and  Accept¬ 
ance  Testing  of  Navigational  Set,  Doppler  AN/ASN-128,  Specification  Number 
AV-IT5800-00 1 A ,  Fort  Monmouth,  NJ,  9  March  1979. 


15-20  KM 


Figure  1 


CP2 


t 


LEG  I-CP2-CPI 
LEG  2-CPI-CP3 
LEG  3-CP3-CPI 
LEG  4- CPI -CP2 


i 

< 
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Graphic  presentation  of  an  ideal  flight  test  course 


5.  In  using  the  AN/ASN  128  to  enter  the  magnetic  variation  of  the 
flight  area,  the  Clark  1866  code  of  CL6  should  be  used  anywhere  within  the 
continental  United  States.  Outside  the  continental  United  States,  other 
applicable  codes  should  be  used.  The  Defense  Mapping  Agency's  Index  of 
Grids,  Datums  and  Spheroids2  should  be  used  to  determine  what  code  is 
appropriate. 

6.  The  installation  and  acceptance  document  suggests  using  the 
average  value  by  flying  three  test  patterns  and  computing  a  mean  value.  If 
good  reliability  is  determined  early  in  the  flight  pattern,  then  one  trip 
around  the  course  should  provide  accurate  data  to  perform  the  necessary 
alignment  calculations. 

7.  The  same  individual  should  determine  the  aircraft's  exact  time 
over  the  check  point  for  recording  and  storage  to  assure  consistency  and 
reliability.  This  will  also  eliminate  another  variable. 

8.  The  alignment  procedure  is  a  precise  and  time  consuming  task. 
The  accurate  alignment  of  the  flux  system  will  yield  greater  system 
performance. 


ALIGNMENT  EQUATIONS  AND  COMPUTATIONS 

This  section  deals  with  equations  and  calculations  necessary  in 
aligning  the  aircraft  flux  system.  The  majority  of  the  data  provided  here 
is  also  contained  in  the  installation  and  acceptance  document.  Some  areas 
have  been  clarified  to  a  greater  extent,  and  some  of  the  acceptance 
document  errors  have  been  corrected.  It  Is  included  to  assist  those  who  do 
not  have  computer  facilities  or  the  background  to  utilize  the  computer 
assisted  methodology  that  will  follow  this  section. 

The  computations  performed  will  use  coordinate  values  supplied  in  the 
installation  and  acceptance  document. 

The  schematic  depicted  in  Figure  2  would  indicate  a  flight  from  CP  1 
to  CP  2  and  return  to  CP  1,  thus  allowing  two  sets  of  data  points  to  be 
generated;  i.e.  actual  (plotted)  and  indicated  (displayed  by  doppler  upon 
overflight  of  plotted  point).  Table  1  shows  the  UTM  actual  and  indicated 
coordinates. 


^Defense  Mapping  Agency  Topographic  Center.  Grid,  datum,  and  spherlod  in¬ 
formation,  Plate  I.  DMA  Stock  No.  DATMX52411C1 .  Code  40220,  Washington,  DC, 
20315,  2  September  76. 
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ACTUAL  COORDINATES 
(16S  EL  1414  5439) 


CP2 


ACTUAL  COORDINATES 
(I6S  EL  OOOO  4025) 


Figure  2.  Initial  test  leg  of  flight  course 


TABLE  1 

UTM  Actual  and  UTM  Indicated  Coordinates 


UTM Actual  Coordinates  UTM  Indicated  Coordinates 


Grid 

Easting 

Northing 

Grid 

Easting 

Northing 

Zone 

Area 

(E) 

(N) 

Zone 

Area 

(E) 

(N) 

16S 

EL 

0000 

4025 

16S 

EL 

0006 

4022 

16S 

EL 

1414 

5439 

16S 

EL 

1438 

5432 

Equations  and  Calculations 


Actual  Distance  Traveled  CPI  — ■  »  CP2 

Dj  -  Distance  traveled 
(m)  -  Meters 

A  ■  Delta;  i.e.,  change  or  difference 


Dj (m)  *  c 


\^a2+b2  XI 0  where 


V(1414)2  +(1414)2 


CP2(E)-1414 

CP2(N)=5439 

a  -  CPl(E)-0000 

b  =  CP1(N)*4025 

AE  *1414 

AN  *1414 

X10 

-  V3998792  X10 


-  20,000m 


*  20km 
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Actual  Bearing  Angle  (Deg) 


meters 

Uhere  y  »  northing  distance  (m) 

meters 

x  »  easting  distance  (m) 

0  -  Tan-1  1 
x 

0  -  Tan-1  1414 
1414 

0  -  45° 


Easting  Error  (meters)  -  £Pt(m) 


i  *  indicated 


tPE(m)  -  (E2i  -  E2)  -  (EU  -  EX ) J  X10 
ePe  (m)  =*  (1438  -  1414)  -  (0006  -  0000)]  X10 

«  (24)  -  (6)]  X10 


180  meters 


Northing  Error 


ePll  (m) 


(N2i  -  N2)  -  (NU  -  NX)J  X10 
(5432  -  5439)  -  (4022  -  4025)]  X10 
(-7)  -  (-3)]  X10 


-40  meters. 
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Using  the  Easting  and  Northing  error  values,  as  well  as  the  actual 
bearing  angle,  we  now  compute  the  Along  Track  Error  (eAT)  and  the  Cross 
Track  Error  (eCT)  in  meters  and  percent. 


Where 

eAT(m)  •  eP||  Co3/3  +  eP[  sin  0 

-  -40  Cos  45°  +  180  sin  45° 

eAT(m)  *  98.980  meters 

eAT(%)  -  lAT(m)  X100  -  98^9 80.  .49% 
DT (m)  20,000 


The  cAT%  should  not  exceed  2.0%  of  DT .  Repair  or  realignment  should 
be  considered  if  error  is  much  more  than  1.0%  of  Dj . 


Cross  Track  Error: 

ECf(m)  -  -ePp|  sin  j3  +  eP{  Cos  0 

-  -(-40)  sin  45°  +  180  Cos  45° 

-  155.54 

eCt (%)  -  c_CT (m)  xiOO 
DT(m) 

-  155»54  X100 

20,000 

eCt(%)  -  .77% 


cC-j%  should  not  exceed  5.0%  of  D-j  .  The  authors  believe  that  the  eCj% 
should  not  exceed  2.5%  of  Dy  in  order  to  consider  the  system  reliable 
enough  for  accurate  navigation. 
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FLIGHT  TEST  DATA  REDUCTION  AND  ERROR  DETERMINATION 


This  section  will  present  the  calculation  methodology  for  reducing  the 
flight  data.  All  the  equations  used  are  taken  from  the  Installation  and 
Acceptance  Testing  Document. 

1.  Calculate  the  azimuth  angle  flown  on  each  leg  by  computing  the 
change  in  northing  AN  and  easting  AE.  In  order  to  compute  this  angle  in  the 
correct  direction,  we  must  use  the  following  values: 


If  AE  >0  then  6-90°  -  Tan-1  AN  +  Magnetic  Variation 

AE 

If  AE  <0  then  6  -  270  -  Tan-1  AN  +  Magnetic  Variation 

AE 

2.  If  using  more  than  one  set  of  values,  group  the  azimuth  angles 
into  north,  south,  east,  and  west  groups. 

3.  Determine  the  error  for  each  leg  by  subtracting  the  calculated 
true  azimuth  from  doppler  read  values* 

i.e.  e  -  0  -  0  true 

4.  Calculate  the  mean  error  for  each  leg;  i.e.  (N,S,E,W): 

c  =»  A  2  £ 

N 

5.  Once  the  mean  is  calculated  for  each  leg,  determine  that  the 
error  (e)  in  each  leg  is  not  greater  than  1°  from  the  calculated  mean  (e). 
The  acceptance  document  recommends  that  the  mean  be  recalculated  leaving 
out  the  leg  with  the  mean  error  difference  of  >1°,  or  reflying  the  course 
to  obtain  a  new  set  of  data  points. 

6.  Calculate  the  index  error  by: 

£  I  -  ?N  +  gS+  £E  +  £W 
4 

7.  Calculate  the  13/S  and  E/V3  compensation  displacement  by: 

£  N/S  - 

2 

and  c 

2 

8.  The  value  calculated  for  the  index  error  (eI)  must  now  be  cor¬ 
rected  from  the  AN/ASN  43.  Accuracy  to  be  obtained  should  be  within  +0.1°. 


9.  The  N/S  and  E/W  compensation  displacement  error  should  be  cor¬ 
rected  to  within  +  0.1°. 

One  will  quickly  realize  that  the  alignment  procedure  is  a  very  time 
consuming  process*  Not  only  must  one  spend  many  hours  in  flight,  but  the 
calculations  required  to  determine  system  error  and  compass  alignment  and 
compensation  values  are  arduous  and  prone  to  error. 

In  the  process  of  installing  and  aligning  two  doppler  systems,  the  US 
Army  Human  Engineering  Laboratory  developed  a  computer  program  that  will  do 
all  calculations  and  print  out  all  information  described  thus  far. 

NOTE:  All  calculations  are  based  on  a  flight  within  one  100,000  meter 
square.  If  the  flight  area  should  Include  more  than  one  100,000  meter 
square,  one  should  take  this  into  account  when  entering  coordinate  values. 
Neither  program  will  recognize  a  change  in  the  100,000  meter  square. 


THEORY  AND  OPERATION  OF  THE  PROGRAM 


The  program  is  written  in  Basic  for  the  Hewlett  Packard  Model  HP-9830A 
programmable  desktop  calculator.  The  calculator  will  calculate  and  print 
pertinent  data  for  calibration  purposes.  A  copy  of  the  Basic  program  is 
Included  in  Appendix  A.  A  sample  of  the  Basic  program  printout  is  contained 
in  Appendix  B. 

It  is  assumed  that  you  are  flying  to  three  points  in  a  manner  similar 
to  the  figure  below  (Figure  3). 

3092,1844  5040,2190 

MAG  VAR  W.10.7  MAG  VAR  W.II.O 

CP3  CPI 


CP2  START  8  END 
4660,0744 
MAG  VAR  W.  10.9 

Figure  3.  Flight  test  course  and  associated  coordinate  and  magnetic 
variations. 

The  easting  and  northing  coordinates  and  the  magnetic  variance  associated 
with  the  three  flight  points  must  be  known  (Figure  3). 
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The  flight  pattern  is  as  follows:  Take  off  from  point  2,  proceed  to 
point  1,  proceed  to  point  3,  return  to  point  1,  and  return  to  point  2  where 
the  flight  ends.  As  the  flight  proceeds  to  the  various  points,  the  doppler 
coordinates  should  be  recorded.  The  coordinate  values  will  probably  be  dif¬ 
ferent  from  the  reference  coordinates.  The  doppler  and  reference  coordi¬ 
nates  should  be  set  up  in  a  table  (Table  2)  for  ease  in  entering  the  needed 
data  into  the  computer  for  processing. 


TABLE  2 

Doppler  and  Reference  Coordinates  Setup 


Plotted  Reference  Recorded  Doppler 

Points  Points 


Point 

Easting 

Northing 

Easting 

Northing 

Magnetic  Var. 
of  End  Point 

2-1 

4860 

0744 

5010 

2124 

11.0 

1-3 

5040 

2190 

3187 

1935 

10.7 

3-1 

3092 

1844 

4930 

2274 

11.0 

1-2 

5040 

2190 

4840 

0833 

10.9 

Table 

2  relates 

directly  to 

the  flight 

path  taken. 

This  information 

keyed  into  the  computer  for  processing  but  in  a  slightly  different  line 
order.  When  the  program  is  executed,  the  program  will  prompt  the  user  four 
times  for  input  data.  A  typical  prompt  is  listed  below: 


Input  PT1  to  PT2  &  Magnetic  Variation 


The  user  would  type  in  the  known  easting  and  northing  coordinates  for  point 
1,  then  the  doppler  easting  and  northing  coordinates,  and  the  magnetic 
variation  associated  with  point  2  (the  point  one  is  flying  to).  The  input 
stream  is  keyed  in  via  the  keyboard  followed  by  an  execute  (carriage  re¬ 
turn).  The  user's  response  to  the  prompt  is  taken  from  Table  2  and  appears 
below: 


5040, 2190, 4840, 0833, 10. 9(ex) 
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The  user  must  take  care  to  input  the  data  correctly  according  to  the 
prompts.  A  sample  of  the  four  computer  prompts  and  the  proper  user  respon¬ 
ses  are  given  in  the  exact  order  as  they  will  appear  below.  The  responses 
are  based  on  sample  data  contained  in  Table  2. 

(Computer  Prompt)  Input  PT1  to  PT2  &  Magnetic  Variation 
(User  Response)  5040, 2190, 4840, 0833, 10. 9(ex) 

(Computer  Prompt)  Input  PT2  to  PT1  &  Magnetic  Variation 
(User  Response)  4860, 0744, 5010, 2124, 11. 0(ex) 

(Computer  Prompt)  Input  PT1  to  PT3  &  Magnetic  Variation 
(User  Response)  5040, 2190, 3187, 1935, 10. 7(ex) 

(Computer  Prompt)  Input  PT3  to  PT1  &  Magnetic  Variation 
(User  Response)  3092, 1844, 4930, 2274, 11. 0(ex) 

A  large  portion  of  the  program  is  made  up  of  output  statements.  The 
program  should  be  suitable  for  a  handheld  programmable  calculator,  such  as 
the  HP-41C,  provided  that  it  has  sufficient  memory.  This  would  be  ideal  for 
use  directly  in  the  field;  however,  no  attempt  has  been  made  to  convert  it. 
Although  the  HP-9830A  is  somewhat  portable,  the  name  is  misleading.  It  is 
more  of  a  minicomputer  than  a  calculator.  It  is  not  battery  operated  and 
runs  on  a  110V  AC  line  making  it  difficult  to  use  In  the  field. 

It  should  be  very  easy  to  convert  the  program  written  for  the  HP-9830A 
calculator  into  standard  Basic.  Most  of  the  statements  do  not  require  modi¬ 
fication.  The  changes  that  are  necessary  should  be  both  minor  and  obvious. 
Since  everyone  does  not  use  Basic,  the  program  was  converted  to  ForTran. 
The  ForTran  program  was  written  for  the  PDP-11  series  computer  from  Digital 
Equipment  Corporation  under  the  RT-11  operating  system.  This  machine  uses  a 
16-bit  word  length.  There  should  be  no  major  problem  In  running  this  pro¬ 
gram  on  any  other  machine  using  the  ForTran  program.  A  minor  formatting 
problem  might  occur  where  the  real  variable  Array  C  is  used  to  output  some 
character  data  on  the  A4  format.  On  machines  with  larger  word  lengths,  this 
A4  format  should  be  changed  to  accomodate  the  larger  character  storing 
capacity. 

The  ForTran  program  directs  all  messages  and  all  inputs  via  logical 
unit  5.  On  the  PDP-11  series  machine,  this  is  the  consol  terminal.  All  out¬ 
put  is  directed  to  logical  unit  6  which  Is  the  system  line  printer  on  the 
PDP-11  series  machines.  If  your  machine  uses  another  logical  unit  numbers 
for  the  terminal  and  line  printer,  change  the  variables  LO  (Logical  Output) 
and  LI  (Logical  Input)  to  what  they  should  be  for  your  machine. 

The  ForTran  program  prompts  the  user  interactively  as  before,  but  the 
input  must  adhere  to  the  input  format  listed  below.  A  copy  of  the  ForTran 
program  and  its  output  appears  in  Appendixes  C  and  D. 


XXXXX , XXXXX , XXXXX , XXXXX , XX . X( CR) 
05040, 02190, 04840, 00833, 10. 9(CR) 


APPENDIX  A 

THE  BASIC  PROGRAM 
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i0Diis?C "  Ihput^W' ,bc  4' 5  3 

13  INPUT  fSlul  hit  I  IS  /I?  i  VPS’  V^I«NCE"» 

20  DISP  "INPUT  PT9  rn  b»  ’ ^  1 ^  ^  1  ^ 1*  W 1 1 S  J 
25  INPUT  8t2»  1J,»?2T?/Ir  5i^’3ftC'  V**I&NCE"» 

38  DISP  “INPUT  PTi  ro3l5i2,33'«2. 43, Pt2,si 

H  SST -?»}>;? «. 

IS  WXT 

«a  rcR  J*3  T0  4 

is  «r  ,jKA  2*  r*"4  ™En  ?3 

H  ns  in 

ill  ®to7|1;!,h'M/,‘ >*«  J* )’Z  1  "N"  "IJ'<-i>'<j‘i).n>«i,'i1,Tf5fN1f!5 

ill  Sl;7i 'U™'*™-™ 7.5J.J7, 

ill « ;?  >•«  o.  5  >88 

ill  «  j;  ,‘i  j:i«  Ji !»)”  J.evllm  J' 8  ,*SIN<«  J'  7  J)  J'w  J’ 2  J>  >*>® 


ill «j;ii  •>•„> 

loo  c5ITE  < 13*180) 

183  fi«  iT?* j4*'  CflLCULfiJI0NS  ‘,2/,3€x,  “T8BLE  i  •• 

190  8$C  66.  fin  la  •ca.-.'r  iNDTPoren  ,.. 

NORTH" 


.  7*  wrviin  I 

183  P«l,74  3* 

}|2  SiT^‘803a"E«ST 

III  foUSTT’8"”"* 

“I  IliljfSe’j;-88  e.^28ii,n"te 
III  !&i?;220"“  m°’ 

225  PRINT 


%  L  'HDLC  l  $/ 

indicrted  magnetic 


fiCTUflL  BEHRING" 


COORDINRTE  DIFFERENCE  DIST«E  RWLE' 


230  n 1 n » 


(Continued) 
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273  FORMAT  /, 36X. "TABLE  2"V 

289  WRITE  <15. 285)  „  „ 

283  FORMAT  13X. "ALONG  TRACK  ERROR" » 19X. "CROSS  TRACK  ERROR 


298  WRITE  < 15.295) 

295  FORMAT  “LEG"»9X.  "<M>". 13X.  "<X>".  17X.  "<M>". 13X. "<*>“»/ 

380  WRITE  <15»305>AC1.103«AC1.  1 1  3.  AC  1  >  12  3»  AC  1 . 13  3 

305  FORMAT  " 1-2" > 7X» F7. 1 . 10X. F6. 2* 1 3X» F7. 1 » 18X. F6. 2 

310  WRITE  <15*315>AC2»10  3»AC2»113.AC2»123»AC2*133 

315  FORMAT  "2-1 “ . 7X. F7. 1  * 10X. F6. 2. 13X. F7. 1 . 10X. F6. 2 

320  WRITE  <15.325)AC3,103.AC3>113.AE3.123>AC3»133 

323  FORMAT  " 1-3" * 7X» F7. 1 » 10X» F6. 2. 13X. F7. 1 . 10X. F6. 2 

338  WRITE  <  15. 335) AC  4.  10  3*  AC  4 » 11  3)  AC 4*  12  3.  AC  4.  133 

333  FORMAT  "3-1 " » 7X» F7. 1 » 10X. F6. 2. I3X. F7. 1 » 10X. F6. 2 

340  FOR  J«1  TO  4 

345  B1«AC J+<-l>t<J+l>.23-AC J.23 

350  B2*AC J+<-l)t<J+l>, 1 3-AC J. 1 3 

355  IF  B2<0  THEM  370 

360  BC J.43»90-ATN<B1/B2)+AC J.53 

363  GOTO  375 

370  BC J.4  3«270-ATN<8VB2)+AC J.5  3 
373  BC  J.  1  3*AC  J.  4  3-AC  J.2  3 
380  BC  J » 2  3=AE  J » 3  3-AC  J  *  1  3 
385  IF  BC  J.  2  3<0  THEN  400 

390  BCJ»3  3*90-ATN<BC  J»1  3/BE  J.2  3XAC  J.5  3 

391  IF  BC  J*  3 3<360  THEN  405 

392  BCJ.3  3-BCJ.3  3-360 
395  GOTO  405 

400  BC  J. 3  3»270-ATN<BC  J* 1 3/BC  J, 2  3>+AC  J*  5  3 

401  IF  BC  J*  3  3<360  THEN  405 
482  BC  J.3  3-BC  J.3  3-360 

405  BC  J.5  3»BC  J.3  3-BE  J.4  3 
410  NEXT  J 
415  PRINT 

420  WRITE  <15.425) 

423  FORMAT  36X. "TABLE  3"./ 

430  A$C 1 » 80 3* "LEG  DIRECTION  DEL  NORTH  DEL  EAST 

433  A$C  50.80  3*"  THETA  TRUE  ERROR" 

440  WRITE  <15.445>A* 

445  FORMAT  B 
450  PRINT 


THETA- 


455  WRITE  <15. 460)BC  1, 1  3. BC  1.23. BC  1.33. BC  1.43. BC  1.53 

460  FORMAT  "1-2".5X. "SOUTH” » 8X> F7. 1 > 5X> F7. 1 »F9. 1 »  ?X»  F6. 1 >  F10.  1 

465  WRITE  <15.4?0>BC2. 1  3.  BC  2. 2  3.  BC  2. 3  3.  BC  2. 4  3.  BC  2.  5  3 

470  FORMAT  "2-l".5X.  "NORTH" . 8X, F7. 1 , 5X» F7. 1 . F9. 1.7X.F6. 1.F10. 1 

475  WRITE  <15.  480  )BC  3. 1  3.  BC  3.23.  BE  3. 33.  BC  3.43.  BC  3.53 

480  FORMAT  "1-3".5X. "EAST" . 9X. F7. 1 . 5X. F7. 1 , F9. 1 . 7X. F6. 1 . F10.  1 

485  WRITE  <15.490)BC4.  1  3.  BC  4. 2  3.  BC  4. 3  3.  BC  4. 4  3.  BC  4, 5  3 

490  FORMAT  "3-1". 5X. “WEST” » 9X» F7. 1.5X.F7. 1»F9. 1.7X.F6, 1.F10, 1.2/ 

495  WRITE  <15.500)<BC2.53-BC l>53)/2 

500  FORMAT  "ERROR-N/S"»5X.F10.2»/ 

505  WRITE  <15.510XBC3.53-BC4.53V2 
510  FORMAT  "ERROR-E/W".5X.F10.2./ 

515  WRITE  <15. 520XBC  1.5  3+BC2.5  3+BC3.5  3+BC4. 53V4 
520  FORMAT  "INDEX  ERROR" . 5X. F10. 2 
525  END 
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CALCULATIONS 
TABLE  1 


ACTUAL 

INDICATED  HAGNETIC 

ACTUAL 

BEARING 

EAST 

NORTH 

LEO 

COORDINATE 

COORDINATE  DIFFERENCE 

DISTANCE 

ANGLE 

ERROR 

ERROR 

1-2 

3040.-2190. 

4840.-  833.  10.9 

14572. 

198.0 

100.0 

1550.0 

2-1 

4840.-  744. 

5010.-2124.  11.0 

14572. 

18.1 

-100.0  - 

•1550.0 

1-3 

3040. -2190. 

3187.-1935.  10.7 

19785. 

270.4 

2050.0 

70.0 

3-1 

3092.-1844. 

4930.-2274.  11.0 

19785. 

90.9 

-2050.0 

-70.0 

TABLE  2 

ALONG 

TRACK  ERROR 

CROSS 

TRACK  ERROR 

LEO 

<H> 

<Z) 

(H) 

(X) 

1-2 

-1505.1 

-10.33 

383.7 

2.43 

2-1 

-1504.4 

-10.32 

384.4 

2.45 

1-3 

-2049.1 

-10.34 

92.3 

0.47 

3-1 

-2048.4 

-10.35 

103.1 

0.52 

TABLE  3 

LEG 

DIRECTION 

DEL  NORTH  DEL  EAST 

THETA  THETA  TRUE 

ERROR 

1-2 

SOUTH 

-1357,0  -200.0 

199.3 

198. 0 

1.3 

2-1 

NORTH 

1380.0  150.0 

17.2 

18.1 

-0.9 

1-3 

EA8T 

-255.0  -1853.0 

272.9 

270.4 

2.2 

3-1 

WEST 

430.0  1838.0 

67.8 

90.9 

-3.1 

ERROR-N/S 

-1.09 

ERROR- 

-E/W 

2.47 

INDEX 

ERROR 

-0.12 
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PAGE  001 


FORTRAN  IV 

0001 

0002 

0003 

0004 

0005 

0006 

1 

0007 

0008 

2 

0009 

0010 

3 

0011 

0012 

0013 

4 

0014 

0015 

0016 

5 

0017 

0018 

0019 

0021 

0022 

0023 

75 

0024 

0025 

0027 

0029 

0031 

0033 

105 

0034 

0035 

115 

0036 

0037 

125 

0038 

0039 

135 

0040 

140 

0041 

0042 

0043 

0044 

0045 

0046 

100 

0047 

0048 

6 

0049 

0050 

7 

VO 2.04 


DIMENSION  A(  4 » 13) «  B  <  4 »  5  ) »C<4) 

DATA  C/' 1-2  ','2-1  ','1-3  ','3-1  '/ 


L  1=5 
LQ=6 

WRITE <  LI ,  1 ) 

FGRMAK1H  ,  '  INPUT  PT1  TO  PT2  &  MAGNETIC  VARIANCE?') 
READ (LI »2)(A(1,J)»J=1»5) 

FORMAT ( 4 (F5. 0, X ) ,F4«1> 

WRITE(LI , 3) 

FORMAT < 1H  ,' INPUT  PT2  TO  PT1  S  MAGNETIC  VARIANCE?') 


READ (LI, 2) <A<2» J) , J=1 ,5) 

WRITE (LI ,4) 

FORMAT ( 1H  ,' INPUT  PT1  TO  PT3  %  MAGNETIC  VARIANCE?') 

READ (LI »2)(A(3»J)»J=1»5) 

WRITE(LI,5) 

FORMAT ( 1H  ,' INPUT  PT3  TO  PT1  g  MAGNETIC  VARIANCE?') 

READ (LI, 2) (A(4,J),J=1»5) 

DO  100  J=1 » 4 

IF ( J, EQ. 2, OR. J, EQ • 4 ) GOTO  75 

A<J,6)=SQRT((A( J,l)-A( J+l»l> >**2+(A( J*2)-A( J+1,2) )**2)*10,0 


A( J+l »6)=A( J»6) 

B1=SQRT( (A( J.2)-A( J+(-l)**( J+l) ,2) )**2) 

B2=SQRT( (A( J»1 >-A( J+(-l )**( J+l ) ,1) )**2) 

IF(A( J+(-l )**( J+l) , 1) ,GT.A( J, 1 >  »AND.A( J+(-l )**( J+l > ,2) .GT. 


*  A ( J , 2 ) ) GOTO  105 

IF(A(J+(-l)**(J+l) ,2) .GT.A(J,2) . AND . A( J+ (-1 >**( J+l ) ,1) ,LT. 

*  A( J, 1 ) )GOTO  115 

IF(A( J+(-l )**( J+l ) ,2) .LT.A( J,2> ,AND.A( J+(-l )**( J+l ) ,1 ) ,LT. 

*  A(J,1))G0T0  125 

IF(A(J+(-l)**(J+l> ,2) .LT«A(J,2) . AND .A(J+(-l)**(J+l) ,1) .GT . 

*  A( J, 1 ) )GOTO  135 

A ( J , 7 ) =ATAN ( B2/B1 ) *57 . 2958+A ( J , 5 ) 

GOTO  140 

A( J, 7 )=ATAN( B2/B1 >  *57 , 2958+A ( J, 5 ) +270 . 0 
GOTO  140 

A(J,7)=ATAN(B2/B1 )*57.2958+A( J, 5) +180.0 
GOTO  140 

A(J,7)=ATAN(B2/B1 ) *57. 2958+A (J, 5) +90.0 

A( J*8)=( (A( J,3)-A( J+(-l >**< J+l ) , 1) )-(A( J+(-l >**(J+1 ) ,3)- 

*  A(J,1)))*10.0 

A( J,9)  =  ( (A( J,4)-A( J+(-l )**( J+l ) ,2) )-(A(J+(-l )**( J+l )  ,4)- 

*  A( J,2> ) )*10.0 

A(J,10)=A( J,9)*C0S(A( J,7)*.017453)+A( J,8)* 

*  SIN(A(J, 7)*. 017453) 

A(J»11)  =  (A(J,10)/A(J*6) ) *  100 . 0 
A(J,12)=A(J, 8) *COS (A(J,7)*.017453)-A(J,9)* 

*  SIN(A(J, 7)*. 017453) 

A(J, 13)*(A( J,12)/A( J,6>  >*100.0 

CONTINUE 

WRITE(L0>6) 

FORMAT ( 1 H  ,34X, 'CALCULATIONS' ,//,36X, 'TABLE  l'/> 

WRITE( LO, 7 ) 

FORMAT ( 1H  ,9X, ' ACTUAL ' , 5X , ' INDICATED ' , 4X, ' MAGNETIC ', 6X , 

*  'ACTUAL' ,2X, 'BEARING' ,4X, 'EAST' ,4X, 'NORTH' ) 


(Continued) 
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FORTRAN  IV 


V02.04 


PAGE  002 


0051 

0052 

0053 

0054 

0055 
0054 
0057 
00 58 
0059 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0074 
0075 
0076 
0077 
0078 
0079 
0081 
0082 
0084 
0085 
0086 
0087 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 


0106 

0107 

0108 

0109 

0110 

0111 

0112 


WRITE<L0»8) 

8  FORMAT ( 1H  . 'LEG' ,X»2<3X, 'COORDINATE' ) ,2X» 'DIFFERENCE' ,4X, 

*  '  DI STANCE '»2X, 'ANGLE' »4X» 'ERR0R'3X, 'ERROR'/) 
URITE(L0,9)C(1),(A<1,J),J*1,9) 

9  FORMAT ( 1H  . A4,F7.0, , F5.0.F7.0, '- ' ,F5. 0,F9. 1 , 7X,F6. 0* 

*  F9.1,F10.1,FB.l) 

WRITE(L0»9)C<2) . <A<2» J) » J*1 ,9) 

WRITE(L0*9)C(3) , <A<3. J), J*1 ,9) 

WRITE<L0»9)C(4)»<A<4»J),J=1»9) 

URITE(LOrlO) 

10  FORMAT ( 1H  ,/,37X, 'TABLE  2',/) 

URITE(LOfll) 

11  FORMAT <1H  ,13X. 'ALONG  TRACK  ERROR' » 19X» 'CROSS  TRACK  ERROR') 
WRITE(L0,12) 

12  FORMAT < 1H  , ' LEG ' » 9X , 2 ( ' (M > ' » 13X , ' < X ) ' * 17X ) / ) 
WRITE<L0,13)C(1),(A(1,J),J=10,13) 

13  FORMAT ( 1H  , A4, 7X,F7 . 1 , 10X,F6. 2, 13X»F7. 1 » 10X,F6. 2) 

WRITE (L0» 13)C<2) , <A(2,J),J=10,13) 

URIT£(L0,13)C(3) *  <A<3, J) , J=10» 13) 

WRITE (LO, 13 ) C(4 ) , < A < 4 > J) » J=10 » 13 ) 

DO  200  J=1 *4 

Bl=A(J+(-l)**< J+1 ) »2)-A< J»2> 

B2=A(J+(-l)**(J+l),l)-A( J»l> 

IF(B2.LT .0.0) GOTO  370 
B<J»4)»90.0-ATAN<B1/B2)*57.2958+A< Ji5) 

GOTO  375 

370  B( J.4)=270.0-ATAN<B1/B2)*57.2958+A( J»5) 

375  B(Jf 1 )=A< J»4)-A( J.2) 

B(J»2)=A(Jf3)-A(J»l) 

IF(B(J»2) .LT .0.0) GOTO  400 

B( J»3)=90»0-ATAN<B< J» 1 )/B< J»2) )*57.295B+A( J»5) 

IF(B(J»3).LT. 360 < 0 ) GOTO  405 
B(J»3)=B(J»3) -360 • 0 
GOTO  405 

400  B(Jf3)=270.0-ATAN(B<J»l)/B(J»2) ) *57 . 295B+A< J r 5 ) 

IF(B<Jf3)«LT. 360 . 0 ) GOTO  405 
B< J»3)*B( Jf3)-360.0 
405  B(J.5)»B(J.3)-B< J.4) 

200  CONTINUE 

URITECL0.14) 

14  FORMAT <1H  ) 

URITE(LO» 15) 

15  FORMAT ( 1H  »36X»'TABLE  3'/) 

WRITE(LO» 16) 

16  FORMAT < 1H  t 'LEG' »3X* 'DIRECTION' »5X» 'DEL  NORTH' »4Xf 'DEL  EAST', 

*  3X, 'THETA', 5X, 'THETA  TRUE' ,3X, 'ERROR'/) 

WRITE(L0,17)C<1), <B(1, J), J=l,5) 

17  FORMAT ( 1H  , A4 , 4X ,' SOUTH' , 8X , F7 . 1 , 5X, F7 . 1 , F9 . 1 , 7X , F6 . 1 , FI 0 . 1 ) 
URITE<L0»18)C<2) , ( B( 2, J) , J=1 , 5) 

18  FORMAT ( 1H  ,A4,4X, 'NORTH' , 8X,F7. 1 ,5X, F7. 1 ,F9 , 1 , 7X, F6 . 1 , F10 . 1 ) 
URITE(L0»19)C<3) , <B(3*J),J=1,5) 

19  FORMAT ( 1H  ,A4,4X, 'EAST' ,9X,F7. 1 ,5X,F7. 1 ,F9. 1 ,7X,F6. 1 ,F10. 1 ) 
WRITE (L0,20)C(4),<B(4,J) ,J=1,5> 

20  FORMAT ( 1H  »A4»4X, 'WEST' ,9X,F7. 1 ,5X,F7. 1 ,F9. 1 ,7X,F6. 1 ,F10. 1 ) 


WRITE<L0,21XB(2,5)-B(1 ,5)  )/2.0 

21  FORMAT ( 1H  , 'ERROR-N/S' ,5X,F10.2/) 

WRITE(L0,22)(B<3,5)-B(4,5) )/2.0 

22  FORMAT ( 1H  , 'ERROR-E/W' ,5X,F10.2/> 

WRITE(L0,23) (B< 1 ,5) +B< 2, 5 )+B< 3,5) +B< 4 ,5) >/4.0 

23  FORMAT ( 1H  , 'INDEX  ERROR' ,5X,F10.2) 

END 
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CALCULATIONS 
TABLE  1 


LEG 

ACTUAL 

COORDINATE 

INDICATED  MAGNETIC 

COORDINATE  DIFFERENCE 

ACTUAL 

DISTANCE 

BEARING 

ANGLE 

EAST 

ERROR 

•M^COh 

(III 

Vt  fO 

5040-  2190 
4860-  744 
5040-  2190 
3092-  1844 

4840-  83 3  10.9 
5010-  2124  11.0 
3187-  1935  10.7 
4930-  2274  11.0 

14572 

14572 

19785 

19785 

198.0 

18.  1 
270.6 
90.9 

100.0 

-100.0 

2050.0 

-2050.8 

TABLE  2 

LEG 

ALONG 

CM) 

TRACK  ERROR 

<7.) 

CROSS 

CM) 

TRACK  ERROR 

a> 

1-2 

2-1 

1-3 

3-1 

-1505. 1 
-1584.4 
-2049. 1 
-2048.6 

-10.33 

-10.32 

-10.36 

-10.35 

383.8 

386.4 

92.5 

183.2 

2,63 

2.65 

0.47 

0.52 

TABLE  3 

LEG 

DIRECTION 

DEL  NORTH  DEL  EAST 

THETA  THETA  TRUE 

ERROR 

*— •  ft)  »— 
1  1  1  1 
•-*  o>~  ro 

SOUTH 

NORTH 

EAST 

WEST 

-1357.0  -208.0 

1380.0  150.0 

-255.0  -1853.0 

430.0  1838,0 

199.3 

17.2 

272.9 

87.8 

198.0 

18.1 

278.6 

90.9 

1.3 
-0.9 
2.2 
-3. 1 

ERROR- 

-N/S 

-1.09 

ERROR- 

-E'W 

2.67 

INDEX 

ERROR 

-0.  12 

NORTH 

ERROR 

1550.0 

1550.0 

70.0 

-70.0 
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